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Abstract 

In this paper we present the analytic form of the heavy flavour coefficient 
functions for polarized deep inelastic lepton-hadron scattering. The expres- 
sions are valid in the kinematical regime Q 2 ^> m 2 where Q 2 and m 2 stand 
for the masses squared of the virtual photon and heavy quark respectively. 
Using these coefficient functions we have computed the next-to-leading or- 
der a s corrections to polarized charm production at HERA collider energies, 
where both the electron and proton beams are polarized. We also give an 
estimate of these corrections at fixed target experiments where the typical 
Q 2 values are much smaller than at HERA. 
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1 Introduction 



Apart from another test of perturbative QCD, deep inelastic electroproduc- 
tion of charm quarks leads to important information about the gluon density 
inside the proton H, 0, @, 0. This is because the dominant production 
mechanism is represented by the (virtual) photon-gluon fusion process || 
which is the only one appearing in the Born approximation. Although in 
higher order of the strong coupling constant a s other subprocesses will also 
contribute it turns out that the above picture remains essentially unaltered. 

Until now almost all attention was paid to unpolarized charm production. 
However in future fixed target || as well as collider |7J , |8j experiments one is 
also interested in charm production in polarized deep inelastic lepton-hadron 
scattering (the case of photoproduction has recently been discussed in || and 
[[UJ). Like in the unpolarized case one is particularly interested in the gluon 
density since it plays an important role in the description of the longitudinal 
spin structure function gi(x,Q 2 ). Here x denotes the Bjorken scaling vari- 
able and Q 2 is the mass squared of the virtual photon exchanged between 
the lepton and the hadron. Contrary to unpolarized charm electroproduc- 
tion, where the cross-section is already calculated up to next-to-leading order 
(NLO), only the Born approximation exists for the polarized case [TT], [O, 
[|Tj|. The calculation of the NLO corrections will be as difficult as that for 
unpolarized charm electroproduction in [Q, which could be only done in a 
semi-analytic way. This is because the cross sections for the parton subpro- 
cesses involve four dimensional phase space integrals where the integrations 
over the azimuthal and polar angles can be performed analytically. The two 
remaining integrations have to be done numerically. After mass factoriza- 
tion the resulting coefficient functions are folded with parton densities so 
that one finally has to do integrations over three variables. Furthermore the 
LO parton cross sections for Fl and F 2 in the unpolarized case are positive 
definite in contrast to those for the spin structure function g\. After mass 
factorization the positive definiteness does not apply to the NLO coefficient 
functions in Fl and F 2 anymore and additional positive and negative parts 
appear for those in g\. This leads to large cancellations in the numerical 
integrations which will particularly complicate the computation of g\. Hence 
the semi-analytic calculation of g\ will be harder than that already carried 
out for Fl and F 2 in [|14[] . Therefore it is important to have an analytic form 
of the heavy quark coefficient functions in some kinematical regime that can 
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serve as a check on the exact 0(a 2 ) calculations which have still to be done 
for the spin structure function g±. Fortunately, as has been shown for Fl 



and F 2 in [IS], one can obtain analytic expressions for the heavy quark co- 
efficient functions in the asymptotic region Q 2 3> m 2 . Here the asymptotic 
formulae for the heavy quark coefficient functions could be inferred from the 
operator matrix elements (OME's) and the light parton coefficient functions 
so that one did not have to resort to cumbersome calculations of loop- and 
phase-space integrals. 

As far as phenomenological applications are concerned, the asymptotic heavy 
quark coefficient functions will be a good approximation at HERA collider 
energies, where both the electron and the proton beams are polarized, be- 
cause there will be events at Q 2 ^> m 2 , where m c is the mass of the charm 
quark. At fixed target energies, where Q 2 is small, this approximation will 
break down but one can partially remedy this by also including threshold 
effects which are due to soft gluon bremsstrahlung. The latter mechanism 
dominates the threshold region of heavy flavour production as is e.g. shown 
in |IB , fTTf for hadron-hadron scattering. By including these threshold ef- 



fects one can obtain a reasonable description of the charm production cross 
section at small Q 2 . We can show this for the unpolarized case since here 
the exact coefficient functions are available [ l4| , flTSfl . Due to the similarity 
between the polarized and the unpolarized cross sections in the threshold re- 
gion one can assume that this approximation will also work for g%{x, Q 2 , m 2 ). 
In this way one can estimate the NLO effects at the much smaller Q 2 values 
characteristic of fixed target experiments. 

The content of the paper can be summarized as follows: in section 2 we 
introduce our notations and give, for polarized Compton scattering, an exact 
analytic expression for the heavy quark coefficient function, which is valid for 
any Q 2 and m 2 . In section 3 we compute the full two-loop spin dependent 
operator matrix elements (OME's) contributing to the spin structure func- 
tion gi{x,Q 2 ). The spin dependent heavy quark coefficient functions will be 
presented in section 4 in the limit Q 2 ^> m 2 . In section 5 we give improved 
expressions for them by including threshold contributions so that they can 
be also used at smaller Q 2 values. Furthermore we make estimates of the 
NLO corrections to polarized charm production at HERA collider as well as 
fixed target energies. The long formulae obtained for the full operator matrix 
elements and the asymptotic heavy quark coefficient functions are presented 
in appendices A and B respectively. 
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2 Heavy flavour production in polarized electron- 
proton scattering 

In this section we present the formulae needed to describe polarized heavy 



flavour electroproduction. Furthermore we summarize the findings in |fL5 
how to derive the asymptotic form of the heavy quark coefficient functions 
from the operator matrix elements and the light parton coefficient functions. 
Heavy flavour production in polarized deep inelastic electron-proton scatter- 
ing proceeds via the following reaction 

e-(4) + Pip) - e-(£ 2 ) + Q( Pl )( Q(p 2 ) ) + 'X'. (2.1) 

Here 'X' represents any final inclusive hadronic state and the momenta of 
the heavy quark (anti-quark), denoted by Q(Q), are given by p\ and p 2 
respectively. The mass of the heavy quark Q(Q) is given by m. Neglecting 
electro-weak radiative corrections the above process is dominated by the 
exchange of one vector boson which carries the momentum q = l\ — i 2 . 
If the virtuality of the exchanged vector boson Q 2 = — q 2 > is not too large 
( Q 2 <C Mf ) reaction (2.1) proceeds via the exchange of one photon only. In 
this case the computation of the cross-section of (2.1) involves the hadronic 
tensor 

W flu (p,q,s) = ^J d 4 ze^ z <p,s\[J^z),J u (0))\p,s>, (2.2) 

where J M stands for the electro-magnetic current and s denotes the spin vector 
of the proton with s 2 = —1 and s.p = 0. The hadronic structure tensor can 
be split into symmetric and antisymmetric parts in the following way 

W, u (p, q, s) = W* u (p, q) + W*(p, q, s) , (2.3) 

W£ v (p, q) = 7^;(SV ~ ^r) F i<( x , Q 2 ) + {p»P» ~ (PnQf + Pv%) 

(p.q) 2 ^ F 2 (x,Q 2 ) ' 

+9i* — — J > ( 2 - 4 ) 

q A > p.q 

W£(p,q,8) = -^W9lAi(^Q 2 ) + - ^)g 2 (x,Q 2 )] , (2.5) 
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where M is the mass of the proton and x stands for the Bjorken scaling 
variable (see below). The structure functions Fl and F 2 already show up in 
unpolarized electron-proton scattering and the heavy flavour contributions 
have been extensively discussed up to next-to-leading order (NLO) in the 
strong coupling constant a s in [14]. If the incoming electron and proton 



are polarized then in addition to Fl and F 2 one also gets the longitudinal 
spin structure function gi and the transverse spin structure function g 2 . The 
latter show up in the polarized electron-proton cross-section 
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dxdydcj) dxdydcj) Q 2 ^ Q 2 ' ^ Q 2 

x[y gi (x,Q 2 )+2g 2 (x,Q 2 )]. (2.7) 

The scaling variables x and y are defined by 

x (0<x<l). y = ^-, (0<y<l). (2.8) 

2p.q pli 

The angle between the spin s of the proton and the momentum £ 2 of the out- 
going electron in (2.1) is denoted by <j>. The lower arrows on a in (2.6), (2.7) 
indicate the polarization of the incoming electron in the direction of its mo- 
mentum l\. The upper arrow on a in (2.6) stands for the polarization of 
the proton which is parallel or anti-parallel to the polarization of the incom- 
ing lepton. The vertical arrows in (2.7) also belong to the proton which is 
now polarized perpendicular (transverse) to the polarization of the lepton in 
either the up or the down direction. In the subsequent part of this paper 
we will limit ourselves to gi(x,Q 2 ) since it contains leading twist two oper- 
ators only whereas g 2 {x,Q 2 ) can also receive twist three contributions. In 
the case of twist two the structure functions can be described by the QCD- 
improved parton model. In this model the heavy flavour contribution to g\ 
can be expressed as convolution integrals over the partonic scaling variable 
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z = Q 2 / (s + Q 2 ) where s is the square of the photon-parton centre-of-mass 
energy (s > 4m 2 ). This yields the following result 



2 v J x z L K z ' q v m 2 fi 2 ' 

^('^flfr, (2.9) 
v z ' y v m 2 ji 2 J > 

where the upper boundary of the integration is given by z max = Q 2 /(4m 2 + 
Q 2 ). The function G(z, fi 2 ) stands for the polarized gluon density. The singlet 
combination of the polarized quark densities is defined by 

nf 

s(^ / i 2 )=^(/ J (^ / i 2 )+7 l (^/i 2 )), (2.io) 
i=i 

where fa and j i stand for the light quark and anti-quark densities of species 
% respectively. The non-singlet combination of the polarized quark densities 
is given by 

rif , n f 

A(z, = 2 [e\ ~ j- E 4) (/*(*, ^) + 7 4 (^ ^)) • (2.H) 
i=i n / fc=i 

In the above expressions the charges of the light quark and the heavy quark 
are denoted by e$ and eg respectively. Furthermore rif stands for the number 
of light quarks and ji denotes the mass factorization scale, which we choose 
to be equal to the renormalization scale. The latter shows up in the running 
coupling constant denoted by a s (fi 2 ). 

Like the parton densities the heavy quark coefficient functions Li (i = q, g) 
can also be divided into singlet and non-singlet parts, which are indicated by 
the superscripts S and NS in (2.9). Furthermore the singlet quark coefficient 
function can be split into 

L* q =L? + L?. (2.12) 
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The above relation originates from the light flavour decomposition of the 
Feynman graphs contributing to the structure function gi(x, Q 2 , m 2 ) in (2.9). 
One class of graphs gives the same contributions to Lg as well as to L^ s 
whereas another class, which we call purely-singlet (PS), only contributes 
to Lg. The latter class is characterized by those diagrams which have a 
gluon in the t-channel and can therefore only contribute to the singlet quark 
coefficient function l£. It turns out that up to order a 2 Lg = L^ s and 
Lg = 0. In the case of the heavy quark coefficient functions Hj there are 
no non-singlet parts and therefore Hg = Hg S . The latter gets contributions 
for the first time in order a 2 , whereas the perturbation series in Hg starts in 
order a s . The distinction between the heavy quark coefficient functions Li 
and Hi can be traced back to the different photon-parton subprocesses from 
which they originate. The functions Lj are attributed to the reactions where 
the virtual photon couples to the light quark, whereas Hi originate from the 
reactions where the virtual photon couples to the heavy quark. Hence Lj 
and Hi in (2.9) are multiplied by ef and e 2 Q respectively. Moreover when the 
reaction where the photon couples to the heavy quark contains a light quark 
in the initial state then it can only proceed via the exchange of a gluon in the 
t-channel between the light and heavy quark (see diagrams 5a, 5b in ||14|| ). 
Therefore this process is purely singlet and there does not exist a non-singlet 
contribution to Hg. 

We will now discuss the various subprocesses which contribute to the 
spin dependent heavy quark coefficient function up to order a 2 . Both the 
coefficient functions Li and Hi can be expanded in a power series of a s / (4tt) 
with coefficients L\ , H^ respectively, where n denotes the order in the 
perturbation series. Until now only the lowest order (LO) coefficient H^ is 
exactly known and it is given by the photon-gluon fusion process 

i*(q) + g(h)^Q(pi) + Q(P2), (2.13) 

which yields the answer [|TTJ, |I2]| , [ |i~3[ | 

H^\z, Q 2 , m 2 ) = T f U(2z - 1) In + 4(3 - Az)s qi ] (2.14) 

y L v 1 — sqi J J 

with 
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where < z < £/(£ + 4) and the colour factor Tf = 1/2 in SU(N). In 
NLO we encounter the following subprocesses. First we have the gluon 
bremsstrahlung subprocess 

7*0?) + 9(h) - g(k 2 ) + Q(p l ) + Q(p 2 ) . (2.16) 

Including the virtual gluon corrections to (2.13) we obtain the contribution 
which is not yet known and for which the asymptotic form in the limit 
Q 2 ^> m 2 will be determined in section 4. In addition to (2.16) we have the 
subprocess 

7*(?) + Q(q)(ki) - q(q)(h) + Q(pi) + Q(p 2 ) , (2.17) 

which can be subdivided into two different production mechanisms. The 
first one is given by the Bethe-Heitler process (see figs. 5a,b in |TJJ]) and 
the second one can be attributed to the Compton reaction (see figs. 5c, d 
in (T^j). In the case of the Bethe-Heitler process the virtual photon couples 
to the heavy quark and therefore this reaction leads to H^ 2 \ Like Hj 2 ^ 
the coefficient function is not known and its asymptotic form will be 
presented in section 4. In the Compton reaction the virtual photon couples 
to the light (anti-) quark and its contribution to L^ s leads to L^ s '( 2 ) = L^ 2 ^ 
(see the discussion below (2.12)). Finally we want to make the comment 
that there are no interference terms between the Bethe-Heitler and Compton 
reactions in (2.17) if one integrates over all final state momenta. 

Like in the unpolarized case the computation of the Compton process is 
rather trivial and we can present an exact form of L^ s '^ 

1 m z zr L ^ 3 1 — z 1 — z y f J)l v z A ' 



+LiL 2 + 2(-DILl + DILI + DIL3 - DIL4)} 

r 8 4 / z \2/ 8 \i 

+ {-3 + — + ((T^) (-16 + 32,-— )}L, 

f 64 112 152 1 / z w512 128 848 2 
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300 2 20 



•7f + 



r l-)} S g 1 ], (2.18) 
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where sqt is defined in (2.15) and C F = (N 2 - l)/(2N) in SU(N). Further 
we have defined 



sq 2 = y 1 - 4| , 

L 1 = ln(l±^) , L 2 = ln(l±^) , L 3 = In f g2 + ^ 

v 1 - sqi ' v 1 - sq 2 1 y sq 2 - sqi 

[1 - z)(l + sqi)\ T .A~sq 2 



DIL, = Li 2 , DIL 2 = Li 2 (±- 

D/L 3 = Li 2 (I^) , D/L 4 = Li 2 (i±^). (2.19) 
v 1 + sq 2 1 v 1 + sq 2 1 



Here Li 2 (x) is the dilogarithm defined in [19 



The derivation of the asymptotic formulae for the heavy quark coefficient 
functions has been given for the spin averaged structure functions Fl and F 2 
in [ 15 1 . The procedure can be immediately carried over to the spin structure 



function g 1 and we only quote the results. 

In the limit Q 2 ^> m 2 the heavy quark coefficient functions Hi can be 
obtained as follows 

with k,£ = q,g and (g> denotes the convolution symbol 

(/ ® 9) (z) = ^ d Zl f 1 dz 2 S(z - z 1 z 2 )f{z 1 )g{z 2 ) . (2.21) 



Notice that for convenience we have suppressed the ^-dependence in (2.20). 
The operator matrix elements (OME's) Am are defined by 



2 

A ^Ct) =<Z\Ok\£>, (2.22) 

where Ok (k = q, g) are the composite quark and gluon operators which 
appear in the operator product expansion (OPE) of the electro-magnetic 
current in (2.2) near the light-cone (see e.g. eq.(2.22) in |[L5||). In (2.20) the 



9 



operators are already renormalized and sandwiched between on-shell light 
quark and gluon states indicated by \£ >. The latter leads to collinear diver- 
gences which, however, are removed via mass factorization so that the A^g 
are finite. Finally the OME's contain one heavy quark loop only and they 
will be calculated up to 0(ot 2 s ) in the next section. The objects in (2.20) 
denote the light quark and gluon coefficient functions which contribute to 
the spin structure function g\. They have been calculated up to second order 
in [ 20(1 . Notice that relation (2.20) also holds for the heavy quark coefficient 
functions of type Li. 

Finally we want to make the remark that, like in the case of the coefficient 
functions, the OME's can also be divided into non-singlet (NS), singlet (S) 
and purely-singlet (PS) parts. In particular we have the relation 

A qq = + , (2.23) 

and the origin of A^ is the same as mentioned below (2.12) for the and 
Hg S . Like in the case of H^ 8 the OME's A^ 8 are determined by the graphs 
where two gluons are exchanged between the heavy quark and light quark so 
that only the singlet channel can contribute. 
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3 Calculation of the two-loop spin operator 
matrix elements 



In this section we present the calculation of the spin dependent OME's in 
(2.22) containing one heavy quark loop up to two-loop order. A similar 
calculation has been carried out for the spin averaged (unpolarized) OME's in 
[ PI . The calculation of the spin dependent OME's proceeds in an analogous 
way. The twist two operators O q and O g appearing in (2.22) are defined by 



O 



Ml"'Mn , 
q,r 



X 



i n S 



fi2 



■ D fln —ip(x) 
2 



+ trace terms . 



(3.1) 



O 



trace terms 



(3.2) 



O 



IMl—IXj, 
9 



X 



i n S 



\ x 



+trace terms . 



(3.3) 



The symbol S in front of the square bracket stands for the symmetrization of 
the Lorentz indices [L\ • • • fi n where n denotes the spin of the operator. The 
above set of operators can be distinguished with respect to the flavour group 
SU(rif) in a non-singlet part (3.1), carrying the flavour group generator A r , 
and the singlet parts (3.2) and (3.3). The quark and the gluon field tensors 
are given by ip(x) and F"(x) respectively with F. 



F^T a where T a 



stands for the generator of the colour group SU(N) (N = 3). The covariant 
derivative is denoted by = + igT a A a (x) where A^x) represents the 
gluon field. The operator vertices corresponding to the operators (3.1) — (3.3) 
can e.g. be found in appendix A of [plj| . The heavy quark coefficient functions 
mentioned in the last section require the calculation of the following OME's 
(see (2.22)). Starting with the photon-gluon fusion reaction (2.13), (2.16) we 
have to compute up to second order 



Q9 V.. 2 



-<g\0 



(n) 



9 > 



(3.4) 



where Q stands for the heavy quark with mass m. The OME's corresponding 
to (3.4) are determined by the Feynman graphs in fig. 1 (first order) and 
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fig. 2 (second order). For the Bethe-Heitler reaction (2.17), which starts in 
second order of a s , we need 

AlY n '(4)=<l\0»\l> ■ (3.5) 

Here the two-loop graphs are depicted in fig. 3. Finally we also want to obtain 
the asymptotic form of the coefficient function for the Compton reaction 
(2.17) presented in (2.18). For this we need to calculate 

^{^)=<^\0)k>, (3.6) 

which is derived from the graphs in fig. 4 containing the heavy quark loop 
(Q) contribution to the gluon self-energy. Since the OME's A$ are S-matrix 
elements they originate from the Fourier transform in momentum space of 
the following connected Green's functions. The connected Green's function 

<0\T(A;(x)O^(0)Al(y)\0> c , (3.7) 

corresponds to equation (3.4). Here Oq'"^ stands for the heavy quark com- 
posite operator which is defined in the same way as in (3.2) except that now 
the Dirac field occurring in Oq'"^ n represents the heavy quark. The heavy 
quark composite operator can be also sandwiched between light quark fields 
denoted by if>i(x) 

<0\T^ l (x)O^"m j (y)\0> c , (3.8) 

which correspond to (3.5). Finally the connected Green's function related to 
(3.6) is given by 

< 0|T(^(x)O^-^(0)^(y)|0 > c , (3.9) 

where again if>i(x) represents the light quarks and 0^ r "^ n is the non-singlet 
light quark composite operator in (3.1). In these connected Green's func- 
tions a and i,j are the colour indices of the gluon field A^ and the quark 
fields ip, ip respectively. Before one obtains the S-matrix elements in (3.4)- 
(3.6) the external gluon and quark propagators appearing in (3.7)-(3.9) have 
to be amputated. Further one has to include the external quark and gluon 
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self-energies given by the graphs s, t in fig. 2. One can simplify the above 
connected Green's functions using the property that the composite operators 
are traceless symmetric tensors under the Lorentz group. Therefore one ob- 
tains many trace terms which are not essential for the determination of the 
OME's ■ These trace terms can be eliminated by multiplying the oper- 
ators by an external source J^...^ = A Ml • • • A Mn with A 2 = 0. Performing 
the Fourier transform into momentum space and sandwiching the connected 
Green's function (3.7) by the external gluon polarization vectors e M (p), e v (p), 
one obtains 

f(p)GQ^(p), (3-10) 

where p stands for the momentum of the external gluon in the graphs of figs. 
1, 2 with p 2 = 0. The tensor Gq^ v has the form 

Gq#» = HT> ^^(A.p)"- 1 ^/ . (3.11) 

Proceeding in the same way for the connected Green's functions in (3.8) and 
(3.9) after multiplication by the Dirac spinors u(p, s), u(p, s) one obtains 

u{p,s)G%u(p,s) (3.12) 

and 

u(p,s)G^\ r u(p,s), (3.13) 

respectively. Here A r denote the generators of the flavour group SU(nf) and 
u(p, s) stands for the Dirac spinor corresponding to the light external quark 
in the graphs of figs. 3, 4 with momentum p (p 2 = 0). The tensors Gq and 
G^ become equal to 

G% = A P Q S / n) (e, % a s yi(A. P r~^ l5 , (3.14) 



G?' = if/' (n) (e,— ^^(A.pr-iAV (3.15) 

Notice that the operators appearing in the connected Green's functions (3.7) — 
(3.9) are still unrenormalized so that the OME's A k £ contain ultraviolet (UV) 
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divergences. Besides the UV-singularities we also encounter collinear (C) 
divergences. They can be attributed to the coupling of massless external 
quarks and gluons in the graphs of figs. 2-4 to internal massless quanta. 
Both types of singularities have to be regularized and we choose the method 
of (i-dimensional regularization. This regularization leads to pole terms of the 
type 1/e* (e = d — 4) in the unrenormalized quantities which have to be 
removed via renormalization and mass factorization. Sometimes it is useful to 
distinguish between UV- and C-divergences and, where appropriate, we will 
indicate them by the notation euv an d ec respectively in the subsequent for- 
mulae. In general d- dimensional regularization is the most suitable method 
to regularize all kinds of singularities appearing in perturbative quantum 
field theories since it preserves the Slavnov- Taylor identities characteristic of 
gauge field theories. However, in the case of spin quantities the 75-matrix 
appears together with the Levi-Civita tensor (see eqs.(3.10)-(3.15)) for which 
one has to find a d- dimensional prescription. For the 75-matrix we will choose 
the prescription of 't Hooft and Veltman |Z| which is equivalent to the one 



given by Breitenlohner and Maison in |23]. Since in our calculations only 



one 75-matrix appears, one can show [24|, [£5 that the replacement of ^75 
in (3.14), (3.15) by 

475 = ^wA'VTV (3.16) 



is equivalent to the prescription given in |22], ||23|| . Although the methods in 



P2j|-[p5| are all consistent they have one drawback, namely that the operator 
0^ r (3.1) gets renormalized in spite of the fact that it is conserved. Fur- 
thermore the 75-prescription in (3.16) also affects 0^ r "^ n (n > 1) (3.1) and 
Oft 1 ' "^"(3.2) so that one has to introduce some additional renormalization 
constants to restore the Ward identities violated by (3.16). Unfortunately 
these constants have been only calculated in the literature E| up to two- 



loop order for spin n = 1 and we need the finite OME's A^ for arbitrary n. 
However, as we will show later on, one can avoid the last procedure as long as 
one combines quantities which are calculated using the same 75-prescription. 
This is what happens in the case of the heavy quark coefficient functions 
Hi as presented in (2.20). If the OME's Am and the light parton coefficient 
functions Ck are computed using the same definition for the 75-matrix and 
the Levi-Civita tensor the prescription dependence will cancel in the product 
on the right-hand side of (2.20) providing us with a unique result for H^. 
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The unrenormalized Aq in (3.11) can be obtained as follows. First one 
replaces the 75-matrix appearing in Oq 1 "' 1 " in (3.7) according to prescription 
(3.16). Because of Lorentz covariance the Levi-Civita tensor in (3.11) emerges 
in a natural way. This tensor can be projected out in (i-dimensions using the 
relation 

W $ ■ «0 = ^- 1(d - 2 ) (d -3) ^< A ^'" G ^ A ^ -< 3 - 17 > 

In this way we obtain only Lorentz scalars in the numerators of the Feynman 
integrals which can be partially cancelled by similar terms appearing in the 
denominators. The traces of the Feynman loops in figs. 1, 2 and the contrac- 
tions over dummy Lorentz indices have been performed with the algebraic 
manipulation program FORM @. The calculat ion of many of the scalar 
integrals has been already done in [15] for the spin averaged analogue of Aq^ 



and we can take over those results except for some additional integrals which 
we have computed for the spin case in (3.17). In the discussion of our results 
we will drop the S in and perform the inverse Mellin transformation 

on the OME's so that they become dependent on the partonic variable z 
in (2.15). This implies that all the products are replaced by convolutions 
according to (2.21). 

The one-loop OME Aq , which only gets a non-zero contribution from 
diagram a in fig. 1, can be cast in the algebraic form 

V 2 / «■ euv 



Here S e is a spherical factor defined by 

5 £ = exp{|[ 7B -ln(47r)]}, (3.19) 

which is characteristic of ci-dimensional regularization and contains the Euler 
constant je- The mass parameter \x originates from the dimensionality of 
the coupling constant g (a s = g 2 /(4ii)) in ^-dimensions and should not be 
confused with the renormalization and mass factorization scales. However, if 
one only subtracts the pole terms like in the MS-scheme, the mass parameter 
/i turns into the afore-mentioned scales. The superscript k in A^' denotes 
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the order in the perturbation series expansion of the OME's which can be 
written as 



4 = E©^- 



(3.20) 



The objects = q, g; k = 0, 1, • • •) which we will need for (3.18) and 

the subsequent expressions denote the spin AP-splitting functions which have 
been calculated up to next-to-leading order in plj , [f27j] . In lowest order the 
renormalization group coefficients in (3.18) become 



p(0) 

a (1) - 



8T,[2*-1] 



-5C(2)P«?. 



(3.21) 



The two-loop OME A^ which is determined by the graphs in fig. 2 is given 
by (Al). As has been shown in section 3 of |U|, it can be expressed into the 
renormalization group coefficients as follows 



1 rl 



H=Q 



6/2 



l + - e ^(2)U 



id) 



(3.22) 



In the above expression, where mass renormalization has already been carried 
out, the pole terms e~ k stand for the UV as well as C-divergences so that 
we have put e uv = e C - The last term in (3.22) can be traced back to the 
graphs s, t in fig. 2. which contain the heavy quark loop contributions to 
the external gluon self-energy. Here one sums over all heavy quarks called H 
starting with H = Q and ending with H = t (top-quark). Notice that Q is 
the lightest heavy quark which is the same as the one produced in the final 
state of process (2.1). Further it also appears in the heavy quark operator 
whose OME's are shown in figs. 1, 2. For instance for charm production in 
(2.1) we have Q = c, m c = m and the sum in (3.22) runs over H = c,b, t. In 
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the case of bottom (b) production we have Q = b, = m and the sum runs 
over H = b,t. The contribution to the beta-function from a heavy quark H 
is given by 



A), 



H 



3 } 



(3.23) 



for all species H which implies that /3o,q = /3q,h- The Riemann zeta-function 
C(2) = 7r 2 /6 originates from the heavy quark contribution to the gluon self- 
energy denoted by U(p 2 ,m 2 ). At p 2 = the unrenormalized expression 
11(0, m 2 ) is proportional to (1 + e 2 C(2)/8)/e. 

The other renormalization group coefficients can be inferred from the 
literature pi] , p7fl-|29|l and they are given by 



A) 



11 r 4 t 
Y C A - -n f T f 



2 ( r l-) + -l-, + ^(l-, ) 



P (o) = 8C 



(r 



+ 1 - 2z 



2/3 5(l-z), 

pW = 4C F T/ [2(1 - 2z) ( - 2 ln 2 (l - 2) + 4 In z ln(l - z) - In 2 z + 4((2) 

+16(1 - z) ln(l - z) - 2(1 - 16z) lnz + 4 + 6z 

+AC A T f [4(1 - 2z) ln 2 (l -z)- 4(1 + 2z) ( In 2 z + 2 In z ln(l + z) 

+2Li 2 (-z)) - 8C(2) + 4(1 + 8z) lnz - 16(1 - z) ln(l - z) 

+4(12 - llz)] , (3.24) 
where the colour factor C A is given by C A = N in SU(N). 

A (2) 

The full analytic expression for the unrenormalized Aq g can be found 
in (A.l). Here the polylogarithmic functions Li n (z) and S njP (z) are defined 
in ||19|| . In order to construct the heavy quark coefficient functions in the 
next section we need to renormalize Aq^ (£ = 0, 1) in eqs. (3.18) and (3.22). 
For the coupling constant renormalization we choose a scheme where the 
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heavy quarks H appearing in the sum of (3.22) decouple from the running 
coupling constant a s (/i 2 ) when /i 2 < m 2 H . This renormalization prescription 
completely removes the last term in (3.22) from the OME Aq g . It also implies 
that only the light flavours indicated by nj appear in the running coupling 
constant. For instance in the case of charm production nj = 3 and the c, b 
and t quark contributions are absent in the running coupling constant. 

Furthermore one has to remove the C-divergences occurring in (3.22) via 
mass factorization. The above procedure has been carried out in section 3 of 



[ |l"5|j in the context of the spin averaged OME and we can simply take over the 



algebraic expressions from that paper. In the MS-scheme the renormalized 



OME's A { q] and A§J become 



1 . m 



2 



41 = -^^ + $1 (3-25) 



+{ - \p$ - «>«a + 5«s ® (p's - 4°')} 

+ag» + 2A,ag' + og> ® (P< 0) - P(?') . (3^26) 



Notice that (3.26) is still dependent on the 75-matrix prescription which en- 
ters the CpTf part (see (A.l)). This is exhibited by the splitting function 
Pjg) (3.24) where the C^T/ part still has to undergo an additional renormal- 
ization in order to become equal to the result in [21], [27]. The prescription 

(2) 

dependence also affects olq'. However as will be shown in the next section 
this will be compensated in the construction of the heavy quark coefficient 
functions when we add the same type of term coming from the massless 
parton coefficient function computed in [20 . 



Finally, the unrenormalized as well as renormalized OME's Aq^ (£ = 1, 2) 
satisfy the relation 

~\zA%l(z) = 0. (3.27) 

PS (2) 

The OME Ag q is determined by the graphs in fig. 3. However it appears 
that only diagram a gives a nonzero contribution. The easiest way to compute 
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this graph is by using the standard Feynman parameterization which provides 
us with (see (3.14)) the unrenormalized OME 



+7(4 P « + °S ® + "IT] • (3-28) 

Like in the case of Ag g (3.22) we did not make any distinction between 
UV- and C-singular pole terms e~ k (eyy = ec). The renormalization group 
coefficients are given by (see also (3.21)) 

4 0) = 4C ^ 2 " A > 

pW) = lQC F T f [-(l + 2) In 2 z - (1 - 3z) Inz + 1 - *] . (3.29) 

The analytic expression for the unrenormalized Aq can be found in (A. 2). 
After removing the UV- and C-divergences the renormalized OME reads (see 
section 3 of [ITS)) 



2 

,2 



+{-k p * (1,+ H>^ , } ln 



2 

m 

2 



+«S' (2) -4>n ( ?- (3-30) 



The above expression is represented in the MS-scheme. Like Aq^ in (3.26), 
Aq^ 2/I still depends on the prescription for the 75 matrix. Since P^'^ al- 



ready agrees with the results in [TZI|], |27| it does not need any finite renormal- 
ization and the dependence on the prescription for the 75 matrix only enters 

PS (2) 

via the non-logarithmic term a>Qq ■ While computing the heavy quark coef- 
ficient function the latter dependence will be cancelled by contributions from 
the massless quark coefficient function C^ s '^ in |[20| . 

NS (2,) 

Finally we call attention to the non-sing let OME A l ™£> derived from the 
graphs in fig. 4. These graphs are also computed by Feynman parameteri- 
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zation and we get the unrenormalized result from (3.15) 



iNS,(2) _ q2 f m Y\ 1 f n p (0)\ 
^99,Q - ^ [-^\ - P0,Q^ qq ) 



-{ 



p 



ns,(i; 

qq,Q 



uv 



}+ 



NS(2) 

l gq,Q 



(3.31) 



Notice that in fig. 4 the heavy flavour loop contribution to the gluon self- 



energy, denoted by II (p 



m 



contains the heavy quark Q with mass m only. 



For the construction of the heavy quark coefficient function L q 2 ^ we do not 
need the contribution of the other heavy quarks (as mentioned below (3.22)) 

a (2) a PS f2l 

which have masses larger than m. Contrary to Aq g and Aq q the above 
equation only contains UV-divergences since the heavy quark Q prevents 

a NS ( 2) 

A qq q from being C-singular provided that the gluon self-energy is renor- 
malized in such a way that ITr(0, m 2 ) = 0. Further in the non-singlet case we 
can completely remove the dependence on the prescription for the 75-matrix 

This has to be so be- 



since after renormalization (^4^„ s ) 



qq yspm 



qq /spm avc 



cause the splitting function P qq is the same for the spin averaged and the 
spin dependent non-singlet operators (see [21], [Pq[ , |27[[ ). Therefore after 



removing the 75 prescription dependence A^q 2 ^ becomes the same as in the 
spin averaged case and we get 



,NS,(1) 



99,' 



= C F T f 
16 1 + z 2 

T i-z 



160 / 1 



176 _ 16 



lnz + 5(1 -z)( 



4 32 



C(2) 



(3.32) 



The analytical expression for A^'q^ can be found in (C.5), (C.6) of 
After renormalization the OME becomes 



4 ns,(2) _ r la P (o)\.2 m2 , r 1 p ns,(i)\i m2 
A 99,Q - i - 4 Aw^m j hi — + I - -i qq Q ) In — 

+«S 2) + 7A, Q C(2)Pi 0) • (3.33) 



Summarizing our above results for the OME's we found that the coef- 
ficients of the double and single pole terms can be inferred from the AP 
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splitting functions and the beta-function which are published in the litera- 
ture. In this way we have a check on the calculations of Aij. The non-pole 
terms Oq g (3.26) a^ q '^ (3.30) and cl^q 2 ^ (3.33) cannot be predicted and are 
calculated in this paper. They are needed to compute the order a 2 s spin- 
dependent heavy quark coefficient functions (2.20) up to the non- logarithmic 
term which will be done in the next section. 
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4 Heavy quark coefficient functions 



In this section we present the heavy quark coefficient functions Hi and Li 
(£ = q,g) defined in (2.9) up to order a 2 in the asymptotic limit Q 2 3> m 2 . 
In [|L5[] we showed that, for the spin averaged case, the asymptotic limits 
for Hi and Li can be expressed into the renormalization group coefficients 
appearing in the algebraic expression for the OME's and the massless 
parton coefficient functions C k . This derivation follows from the mass fac- 
torization theorem, which is also applicable to the spin structure function 
gi(x, Q 2 ) in (2.9). The theorem states that the mass dependent terms of the 
type In 1 (m 2 /fi 2 ) In- 7 (Q 2 /m 2 ) occurring in the asymptotic expressions for Hi, 
Li can be factorized out in the following way 

L i&.£)=M%) (4.2) 

where Hi, Li {t = q, g) denote the spin dependent heavy quark coefficient 
functions in the limit Q 2 ^> m 2 and r =S, NS. Up to order a 2 the finite 
spin dependent OME's A ki (k, £ — q,g) have been calculated in section 3 
and are given in Appendix A. The spin dependent light parton coefficient 
functions have been also calculated up to order a 2 and they can be found 



in [|20] . Notice that A ke as well as C k have been computed in the MS-scheme 
and in the case of (4.1) both depend on the prescription for the 75-matrix. 
In the products appearing on the right hand sides of (4.1), (4.2) the scheme 
dependence is only partially cancelled, which is revealed by the dependence of 
Hi, Li on the factorization scale /i 2 . The latter originates from the coupling 
of a light parton (gluon or quark) to an internal light parton which appears in 
subprocesses (2.16), (2.17). Fortunately the 75-matrix prescription cancels in 
the product on the right hand side of (4.1) as we will discuss below. Since the 
algebraic structure of the spin dependent heavy quark coefficient functions 
corresponding to gi(x, Q 2 ) in (2.9) is the same as derived for the spin averaged 
structure function ^(x, Q 2 ) we can simply take over the formulae in [fLUR . For 
the LO photon-gluon fusion process (2.13) we obtain 

iA =I P (o) ln ^ 

-m? ' jj? ' 2 qg m 



Q 2 _rn\_l {0) (1) (1) 

2 ' ,,2 J - 9 19 111 ™2 + a Qg + C 9 ' V 1 -*) 



22 



where and Oq appear in Aq g (3.25) and is the non-log term in the 
lowest order coefficient function defined by (see p0| ) 




Notice that like Am in (3.20) the coefficient functions are expanded as 

(4-5) 

The explicit expression for (4.3) can be found in (5.1). In NLO the spin 
dependent heavy quark coefficient function of the photon-gluon fusion process 
(2.13), (2.16) becomes 

+{\p$ ~ + \p® « «g + \p!S » 4" + 5^ ® 4"} 

+ { " A(c<" + «g) + ±P< a °> 9 (c« + „«)} In J 

+cf + ag» + 2,3„og + cW ® ag» + P<J> 8 ag - P« 8 ag> , 

(4.6) 

where the coefficients P^, P}j \ QqL ®Qg show up in AqI (3.26). The co- 
efficients c^p (£ = 1, 2) show up in the order a 2 s contribution to the gluon 
coefficient function calculated in PD| . The latter can be written as 

+ { \p$ - A4" + 5^ « 4" + ^S> » 4" } i -J 

+4 2> • (4-7) 

The explicit expression for (4.6) is given in (5.2). The latter can be split 
into CVT/ and CaT/ parts (see (B.2)). The C \Tf part still depends on the 
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mass factorization scale /i 2 and is therefore scheme dependent (in our case 
the MS-scheme). 

The same scale dependence was also found for the exact expression of 
(4.6) in |]14j| . It can be attributed to the coupling constant renormalization 



represented by the lowest order coefficient /3 in the beta-function and to 
mass factorization which is revealed by the lowest order splitting function 
Pfjg'- The CpTf part of is scheme independent which implies that it 
is obtained without performing renormalization and mass factorization on 
the original parton cross section of the photon-gluon fusion process (2.13), 
(2.16). Therefore the latter did not contain UV and C-divergences and it 



can be computed in four dimensions. From section 3 and |2(J we infer that 
the prescription dependence of the 75-matrix only enters in the CfTj parts 
of A { q] (3.26) and (4.7). Hence it has to cancel in the same part of 
because the latter can be computed in four dimensions, where the 75- 
matrix is unique. Further we made an interesting observation for the gluonic 



coefficient functions H g and C g (see |2(|) i.e, 



f 1 dz Hf\z, ?t) = , /' dz Cf{z, %) = 0, (4.8) 
Jo wr Jo jj, 



with t = 1,2. This property was already mentioned for Aq g (see (3.27)). 

The asymptotic expression of the order a 2 s spin dependent heavy quark 
coefficient function corresponding to the Bethe-Heitler process in (2.17) is 
given by 

$) = {^;> » } * $ + {\?:, U1) + H 0> s 111 § 

+{\P!S> « P?A in $ ' » J + {\P?1 » (<£> + «g>} 1» 5 
+< S ' (2) +a?/ 2 »-P(»»®ag». (4.9) 

The coefficients -P^ s,< - 1 ' ) , -P/^ , a^''' 2 '' and a!q g show up in A^^ (3.30) and 
Cq S '^ appears in the pure singlet quark coefficient function calculated in 



c f' ,2, (f) = {H° ,SP »} ln2 ^ 
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+{^ s ' (1) + 5^ » 4"} h J + < s ' <2) • (4-10) 

The explicit expression for (4.9) can be found in (B.3). Notice that is 
still scheme dependent (in our case the MS scheme), which is indicated by the 
mass factorization scale /i 2 . This factorization scheme dependence shows up 
via the splitting function Pf$- The 75-matrix prescription dependence enters 

via the nonlogarithmic terms in A^ q '^ (3.30) as well as in C q s '^ (4.10) and 
therefore cancels in the sum of both which equals . 

Finally we turn our attention to the Compton process in (2.17) which 
provides us with the spin dependent coefficient function L^ s . As can be 
inferred from (2.18) it is scheme independent since there is no dependence 
on fi 2 . Furthermore the prescription dependence for the 7 5 -matrix could 
be removed from A^ q S '^ as well as from the order a 2 non-singlet coefficient 



function . The latter reads (see (4.27) in fTSfl) 



r NS,(2)/Q m \_( 1 o P (o)l ln 2<5 

4-I 1 /! P(°)\1ti 2 ^ + 1 P NS 'Wln ^ 



H 2 2 fi 2 
A Q c«lng + C S^. (4.11) 



In the above expression we have chosen the same renormalization for the 
heavy quark loop (Q) contribution to the gluon self-energy as the one ap- 

NS (2) 

pearing in the OME A 'X (see below (3.31)). Finally the heavy quark 
coefficient function is given by 



+{kS (1, -^ I, }'«S+« 2, +« 



+ jA>,«<(2)P<? , (4.12) 



where the eoefflcients P, N , s «f , P») , /3„, e and a v Jf show up in < s f (3.33). 
The coefficient can be found in [20] . The explicit expression for (4.12) 
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can be found in (B.4). One can check that in the limit Q 2 ^> m 2 the exact 
expression for the Compton process in (2.18) becomes equal to the asymptotic 
formula in (4.12) or in (B.4). 
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5 Results 



In this section we compute the heavy charm component of the spin structure 
function gi(x, Q 2 ) and compare it with the light parton contributions. As 
has been mentioned before (see section 2) the order a s spin dependent heavy 
quark coefficient function (2.14) is exactly known whereas the order a 2 
coefficient functions computed in section 4 are only valid when Q 2 3> m 2 . 
Therefore the latter can only be used at those Q 2 -values which are charac- 
teristic for the polarized electron and proton beam facility at HERA J7J, [|]. 
If we want to use the asymptotic expressions at the lower Q 2 -values, which 
are typical for the SMC-experiment M, we have to improve them. To this 
purpose we will construct below some improved coefficient functions which 
can be also used at smaller Q 2 -values. To check the validity of this approx- 
imation we carry out the same procedure for the spin averaged structure 
function F 2 (x, Q 2 ,m 2 ). The latter is expressed in parton densities and co- 
efficient functions in the same way as done for gi(x,Q 2 ,m 2 ) in (2.9). The 
heavy quark coefficient functions needed for both gi and F 2 start in the same 
order of a s and they have the same asymptotic behaviour exhibited by the 
large logarithmic terms ln*(m 2 //z 2 ) \n J (Q 2 / m 2 ) . Further the spin averaged as 
well as the spin dependent Hp* show the same threshold behaviour which is 



determined by soft gluon bremsstrahlung | 16[] , ||17|| . However contrary to g\ 
the exact heavy quark coefficient functions are known for F 2 and are pub- 
lished in ||14|| . Therefore we can test the quality of the approximation for F 2 



and assume that it also works for g± using the same values for Q 2 . Since we 
already have the exact form of the coefficient functions L^ s '^ corresponding 
to the Compton process in (2.17) there is no need for an approximation here. 
The same holds for the Born approximation given by in (2.14). 

In the case of the bremsstrahlung reaction (2.16) and the Bethe-Heitler 
process (2.17), which lead to the heavy quark coefficient functions and 
respectively, we propose the following approximation 

+ (-)VM)^ s ( S ,m 2 ) (5.1) 
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with 



SthresO,?™ 2 ) = C A 4 In 2 (l 



4m z 



201n(l 



Am" 



+4 In — In ( 1 - 

ii, 2 \ 



C 



Am 2 



2tt 2 



Am 2 / s 



(5.2) 



and 



Am 2 \ 1/2 



(5.3) 



where £ and 2 are defined in (2.15) and s is the virtual photon-parton CM 
energy. Expressed in z and £ the latter becomes equal to 



1 - z 



(5.4) 



The spin dependent heavy quark coefficient functions on the right hand side 
of (5.1), (5.3) are given in appendix B and they are strictly valid for £ = 
Q 2 /m 2 3> 1. To improve their behaviour near threshold (s = Am 2 ) we 
multiply them by the factor (1 — Am 2 /s) 1 ^ 2 . Further we add to (5.1) a 
term which is obtained by multiplying the exact Born coefficient function 
Hy> (2.14) by the factor S^res (5.2). The logarithmic terms in this factor 
originate from soft gluon bremsstrahlung which is the dominant production 

(2.16) 



mechanism near the threshold in process (2.16) (see ||16|| , Jl7||). The last 
term in (5.2) represents the Coulomb singularity which can be attributed to 
the loop-graphs where one gluon is exchanged between the massive quark 
antiquark pair. The factor S t hres is universal (see |17| ) and is the same for 
F 2 and g%. It has been computed for F 2 in eq. (5.7) of fI3fl . To suppress 



unwanted effects at larger s we have removed the factor 8 in the argument 
of the logarithms in eq. (5.7) of |I4J and multiplied S^res with (4m 2 /s) 1//4 so 
that the second part in (5.1) vanishes when s 3> m 2 . 

To test the approximation made above we apply it first to F 2 (x,Q 2 ,m 2 ) 
for which the exact [TJI] as well as the asymptotic coefficient functions H^i 
(I = q,g), valid for Q 2 3> m 2 , (see appendix D in fl5|l ) are known. To that 
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purpose we plot the ratio 

«('.^ = F£ff£) (5 ' 5) 

in NLO and examine below for which £ (or Q 2 ) this approximation breaks 
down. Here F 2 is given by the same formula as in eq. (2.9) where now 
the heavy quark coefficient functions stand for the spin averaged ones. In 
eq. (5.4) the superscripts 'exact' and 'approx' indicate whether the exact 
heavy quark coefficients in Jl4| are used or the spin averaged analogues of the 
approximations in (5.1), (5.3). Notice that for the Born coefficient functions 
the exact formula has been taken. 

For our plots of R in (5.5) we have chosen the parton density set GRV in 
the MS-scheme f3(|. Further we limit ourselves to charm production (m c = 
m = 1.5 GeV/c 2 ) which implies that the number of light flavours rif in (2.9) 
has to be taken to be three in the coefficient functions and the running cou- 
pling constant a s (fi 2 ) (A^ s = 248 MeV). The factorization/renormalization 
scale /i 2 is set equal to Q 2 . We have studied R(x,Q 2 ,m 2 ) in the range 
1(T 4 < x < 1 and 10 < Q 2 < 10 4 (GeV/c) 2 . 

In fig. 5 we have plotted R (5.5) function of Q 2 at x = 0.1, 0.01, 10" 3 
and 10~ 4 . From this figure one infers that for Q 2 > Q^m > 20 (GeV/c) 2 R is 
very close to one (actually 0.9 < R < 1.1) which means that above this value 
^exact an( j ^approx co j nc j ( j e _ E ven for Q 2 > 10 (GeV/c) 2 the approximation 

is rather good if one bears in mind that the statistics of deep inelastic charm 
production is quite low. For x > 0.1 the approximation gets worse which is 
revealed by fig. 6. Here R > 1.2 which happens for x = 0.2 when Q 2 = 10 
(GeV/c) 2 or x > 0.5 when Q 2 = 100 (GeV/c) 2 . Therefore we can conclude 
that the approximation breaks down at large x and small Q 2 -values. In the 
case of the HERA collider this is not bad because the large x-regime is not 
accessible. However for fixed target experiments where x is large and Q 2 is 
small our predictions of the NLO corrections have to be considered as an 
order-of-magnitude estimate only. 

Next we study the validity of our approximation for the charm compo- 
nent of the spin structure function denoted by gi(x, Q 2 , m 2 ) (2.9). Here 
we choose the leading log (LL) parametrization and the next-to-leading log 
(NLL) parametrization in the MS-scheme for the spin parton densities in 
[j3l"H . Here one has two sets of parton densities which are obtained in the 
standard scenario and the valence scenario. Our subsequent plots are made 
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in the standard scenario although we also studied the valence one. It turns 
out that the differences between both scenarios are irrelevant for charm pro- 
duction so that we do not present separate figures for the valence scenario. 
Further the number of flavours and the running coupling constant are the 
same as those taken for F 2 (x, Q 2 , m 2 ) above. Since the exact coefficient func- 
tions f/^ 2 )' cxact — are unknown a comparison between the exact and 
approximate spin structure function g\ can be only made on the Born level. 
Furthermore g\ is not a positive definite quantity either in the exact or ap- 
proximate formulae. This implies that the numerator and the denominator 
in (5.5) can vanish so that it makes no sense to plot R in the polarized case. 
Therefore we have to study the approximation on the level of the structure 
function itself. 

In figs. 7a, 7b and 7c we have plotted gf""* and g approx on the Born level 
for Q 2 = 10,50 and 100 (GeV/c) 2 respectively. Here # approx is obtained from 
the asymptotic coefficient function in (B.l) by multiplying the latter with 
(1 - 4m 2 /s) 1/2 . At Q 2 = 10 (GeV/c) 2 (fig. 7a) the deviation between the 
exact and approximated result is of the order of 25% for x < 0.02. For 
larger x-values the approximation becomes much better. When Q 2 increases 
(see figs. 7b, 7c) g^ xact (Born) and g approx (Born) almost coincide over the 
whole x-region. In NLO we have to take the approximate order a s correction 
to the charm structure function g±(x, Q 2 , m 2 ) because not all exact order a 2 
heavy quark coefficient functions are known. Denoting this approximation by 
<7 approx (NLO) we assume that its validity holds for the same x and Q 2 -values 
as those observed for F 2 approx above i.e. x < 0.1 and Q 2 > 10 (GeV/c) 2 . For 
i > 0.1 we expect that the exact NLO charm structure function will be very 
small which is already indicated by the Born contribution in figs. 7a, 7b and 
7c. The expression for g approx (NLO) is obtained by adding to gJ xact (Born) 
the order a s corrections. The latter originate from the exact expression for 
Lf'( 2 ) (2.18) and H^ appiox (I = q,g) given in (5.1), (5.3). The results for 
^approx ^ ^ o are presented in figs. 8a, 8b and 8c for Q 2 = 10, 50 and 100 
(GeV/c) 2 respectively, where they are compared with the Born contribution. 
From these figures we infer that the bulk of the correction occurs in the region 
0.01 < x < 0.1 and amounts to almost 100% ( # appr °7#f act (Born) ~ 2). 

In figs. 9a, 9b and 9c we have plotted in NLO gt Pprox (x,Q 2 ,m 2 ) and 
gi ght (x, Q 2 ) where the latter structure function is due to light partons (u,d,s 
and g) only. Below x = 2 x 10~ 3 the charm as well as the light parton 
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contribution can become negative so that in this region we have taken the 
absolute values of g\. Further we observe that the charm component of the 
total structure function c/f ght + gf w °* is small. At Q 2 = 10 (GeV/c) 2 (fig. 9a) 
it becomes maximal at x = 10~ 3 where it amounts to 14% of the light parton 
contribution. At larger Q 2 -values i. e . Q 2 = 50, 100 (GeV/c) 2 (see figs. 9b, 9c) 
the maximum occurs at x — 0.01 where it is only 4%. From these results 
we conclude that the charm component of the structure function gi is much 
smaller than the one discovered for in |14]] . In the latter case it becomes as 
large as 40% in the small x-region (x = 10~ 4 ). The above predictions made 
for gi(x,Q 2 ,m 2 ) up to NLO can be considered as reliable because the bulk 
of the corrections are in the region x < 0.1 at Q 2 - values for which ^ pprox 
and g^ pprox (Born) are close to their exact values. This region is accessible 
to fixed target and HERA experiments where in the latter case both the 
electron beam and proton beam are polarized. 

Summarizing our work we have computed the order a 2 contributions to 
the heavy quark coefficient functions corresponding to the spin structure 
function gi(x,Q 2 ,m 2 ). For the Compton process (2.17) we were able to 
calculate the exact coefficient function L^ s '( 2 \ In the case of the photon- 
gluon fusion process (2.16) and the Bethe-Heitler process (2.17) we could 
only obtain an analytic expression of ifj (I = q, g) in the asymptotic regime 
Q 2 > m 2 . The last expressions will serve as a check on the exact coefficient 
functions which can be computed in a semi-analytic way similar to the pro- 



cedure outlined in |H |. To estimate the NLO corrections to gi(x,Q , m ) 



we modified the asymptotic coefficient functions according to (5.1), (5.3) in 
order to mimic the exact form. This approximation was tested for F2 and 
gi (Born) and lead to reasonable results as long as Q 2 > 10 (GeV/c) 2 and 
x < 0.1. Using this approximation we found that the order a s corections to 
gi(x,Q 2 ,m 2 ) are large but the charm component of the total spin structure 
function given by (?i ght + g^ pprox is small. 
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Appendix A 



Here we present the unrenormalized spin dependent operator matrix ele- 
ments A\j whose general structure was expressed in renormalization group 
coefficients in section 3. After having carried out mass renormalization the 
two-loop OME is given by the following expression (see also (3.22)) 



aW(— A - 9 



^{C F 7>[(16-32z)[lnz-21n(l - z)\ +24 
+C A T f [(32 - 6Az) ln(l - z) - 64(1 + z)\nz- 192(1 - z)}} 
+^{C F T f [(4 - 8z) [2 ln 2 (l - z) + In 2 z - 4 In z ln(l - z) - 4£(2)] 
-32(1 - z) ln(l - z) + (4 - 64z) In z - 8 - 12z 
+C A T f [(8 + 16z)[2U 2 (-z) + 21nzln(l + z) + In 2 z] 
-(8 - 16*) ln 2 (l - z) + 16C(2) + 32(1 - 2) ln(l - z) 
-(8 + 64z) In 2 - 96 + 88z] } + a§ g {z) 



H=Q 



[I 



x 



(» + 7« 2 >) 



(A.l) 



with 



ag>(z) = C F 7>{(-1 + 2z)[8C(3) + 8((2) ln(l - z) + - ln 3 (l - z) 

-81n(l -*)Li 2 (l -2) +4C(2)ln^-41nzln 2 (l -2) 

+ - In 3 2 - 81nzLi 2 (l - 2) + 8Li 3 (l - z) - 24S 1)2 (1 - 2)] 
3 

-(116 - 48z - 16z 2 )Li 2 (l - z) + (50 - 32z - 8z 2 )C(2) 
-(72 - 16z - 8z 2 ) In z ln(l - 2) + (12 - Sz - Az 2 ) ln 2 (l - z) 
-(5 - 8z - 4z 2 ) In 2 z - (64 - 60z) ln(l - z) 
-(16 + 50z) lnz- 22 + 46z} 

+C A T f {(-l + 2z)[-8C(2) ln(l - z) - \ ln 3 (l - z) 
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z 



-8C(2) ln(l + z) - 161n(l + z)Li 2 (-z) - 8 ln^ln 2 (l + 2) 

+4 In 2 2 ln(l + z) + 8 In *Li 2 (-z) - 8Li 3 (-*) - 16Si )2 (-z)" 

+ 16(1 + ^)[4Si )2 (l - z) + 2 In*Li 2 (l - 2) - 3((2) Ins + Li 2 (-z) 

+ lnzln(l + z)} - 16(1 - z)C(3) + (100 - 112* - 8s 2 )Li 2 (l - z) 

-(132 - 144z - 4^ 2 )C(2) - 4^(4 + z)\nz ln(l - 2) 

-(10 - 8z - 2z 2 ) ln 2 (l - 2) - (6 + 2s 2 ) In 2 2 

+4 ln(l - 2) - (56 + UBz) \nz- 204 + 212,2 } . (A.2) 



The unrenormalized OME corresponding to fig. 3 (see (3.28)) is given by 



= S^)'CrT,{U - 32(1 + ,)ln, - 80(1 - ,)) 

Lb Lt c 



aQ g A2) (z) = (l + z)[32S li2 (l-z) + 161nzLi 2 (l-z) -24C(2)lnz 




(A.3) 



with 




(A.4) 
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Appendix B 



In this appendix we present the spin dependent heavy quark coefficient func- 
tions iff and Lf^ (i = q,g) in the asymptotic limit Q 2 ^> m 2 . Starting 
with the lowest order photon-gluon fusion process (2.13) the heavy quark 
coefficient function reads (see (4.3)) 



^ (*' T/h4(1 " 2Z) ( M S + 111(1 " Z) ~ M Z ) 



+4(3 - Az)] 



(B.1) 



In next-to-leading order the coefficient function corresponding to the vir- 
tual corrections to the Born reaction (2.13) and the bremstrahlung process 
(2.16) are given by (see (4.6)) 

tff(^'^r) = [C F T f {(8z-4)(2Hl-z)-\nz) + 6} 



Q 2 

2 



m 



+C A T / {-8(l - 2z) ln(l -z) + 16(1 + z) \nz + 48(1 - z)}\ In 2 

+ [C F T f {-B(l -2z)[Li 2 (l -z) -3mzln(l-z) 
+21n 2 (l-^) +ln 2 z-4C(2)] 

+4(17 - 20z) ln(l - z) - 16(3 - 2z) In z - 4(17 - IZz)} 
+C A T){-16(1 + 2z)[U 2 {-z) + lnzln(l + z)\ 
+32(1 + z)U 2 (l - z) + 48 \nz ln(l - z) - 8(1 - 2z) ln 2 (l - 2) 
-8(3 + Az) In 2 z - 32zC(2) + 16(7 - 82) ln(l - z) 

Q 2 



-24(5 - Az) In z - 8(20 - 2lz)} 



In 



r O 2 

+C A T f {-16(1 - 2z) ln(l - z) + 32(1 + 2) In z + 96(1 - z)} In 

+32(1 + z)Li 2 (l - 2) + 48 lnzm(l - 2) - 16(1 - 2z) ln 2 (l - 2) 

2) 



-16(1 + z) In 2 2 + 16(9 - 10*) ln(l 
-32(4 -z)\nz + 16(1 - 2z)C(2) - 256(1 - z) 



2 
/i 2 



(1 - 2z)[24Li 3 (l - z) - 81n(l - z)Li 2 (l - 2) 
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— 32C(2) \nz - 81n 3 (l - 2) + 20hi2ln 2 (l - z) 

o 

-161n 2 2ln(l-2) + -ln 3 z\ 
o 

-16(l + 2) 2 [4S lj2 (-2) +41n(l + z)Li 2 (-2) 
+2 In z ln 2 (l + z) - In 2 2 ln(l + 2) + 2£(2) ln(l + z)\ 
-(32 - 1922 + 322 2 )Li 3 (-2) - (96 - 162 2 )Li 2 (l - z) 
+32(1 - 2) 2 [S li2 (l - z) + lnzLi 2 (-z)] 

+ (^ + 642 + ^2 2 ) [Li 2 (-*) + In z ln(l + 2)] 

oZ o 

+ (66 - 8O2 - 42 2 ) ln 2 (l - z) - 32z 2 ((2) ln(l - z) 

92 

-(188 - 1642 - I62 2 ) ln(l - z) + (36 - 8z z 2 ) In 2 z 

-(160 - 112* - 82 2 ) lnzln(l - z) + ^(320 - 4242 - 482 2 ) In 2 
-(48 - 224,2 - 32^ 2 )C(3) - ^(192 - 3362 - 1842 2 )C(2) 
+^(304- 244^) 



+CU7>[l6(l + 2z)(Uz( L — -) - Li 3 ( - y 

X 1 ,2- X 



1-z 



+ Z' 

- ln(l - z)U 2 (-z) - lnzln(l - 2) ln(l + 2)) 

+8(1 + 2z + 22 2 ) [2Si, 2 (-z) + In z ln 2 (l + 2) 

+2 ln(l + z)Li 2 (-z) + C(2) ln(l + 2)] 

+(72 + 8O2 - 16^ 2 )Si i2 (l - z) - 16(2 + *)[2Li 3 (l - z) 

+ In 2 2 ln(l - z)} - (12 + 242 - 82 2 ) [2Li 3 (-z) 

-ln 2 zln(l + 2) - 2\n zU 2 (-z)\ 

- V f" + 3 + 92 + n ^ 2 ) [Li 2 (-^) + In z ln(l + 2)] 

3 K z ' 

+24 In 2 In 2 (1 - z) + 161nzLi 2 (l - 2) 

+32(1 + z) ln(l - z)Li 2 (l - z) + (76 - 1122 - 82 2 )Li 2 (l - 2) 
+ (38 - 482 + 2z 2 ) ln 2 (l - z) + (24 - 8O2 + 162 2 )C(2) ln(l - z) 

+jj(3 + 4^) In 3 z + (84 - 32z + y2 2 ) In 2 2 

-(136 - 1122 + 42 2 ) hi2ln(l - 2) - (80 + 322)C(2) \nz 
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+-(776 - 652z + 24z 2 ) lnz - (172 - 212z + 8z 2 ) ln(l - z) 

-(28 + 242 + 16z 2 )C(3) - (228 - 224z + ^z 2 )((2) 

+ i(808 - 832z)l . (B.2) 
o 

The coefficient function corresponding to the Bethe-Heitler process (2.17) 
reads (see (4.9)) 



H ? (*• $ ) = C - T l [{»( 1 + *> "» * + 20( 1 - *) } In 2 g 

o 2 

+{ [16(1 + 2) In ^ + 40(1 - z)] In 

+8(1 + z) [2Li 2 (l - 2) + 2 In z ln(l - z) - In 2 z] + 40(1 - z) 

2 

777, 

xln(l -z) - (56 -8z) lnz - 96(1 - z) ) In — 
+{16(1 + z)[Li 2 (l - z) + lnzln(l - z) - In 2 z] 

O 2 

+40(1 -z) ln(l - z) - 32(2 - z) lnz - 88(1 - z)}ln ^~ 

+ (1 + z)(32Si, 2 (l - z) - 16Li 3 (l - z) + 81nzln 2 (l - z) 

-161n 2 zln(l - z) + 161n(l - z)Li 2 (l - z) - 32C(2)lnz 
16 

+y In 3 z) + 16(1 - 3z)Li 2 (l - z) 

-(p + 32 + 32z + y z 2 ) [Li 2 (-z) + In z ln(l + z)\ 

-(112 - 80z + y^ 2 )C(2) - 32(2 - z) lnzln(l - z) 

/ 592\ 
+ (1 - z) (20 ln 2 (l - z) - 88 ln(l - z) + —J 

+ ( 56 + l A z ^ i n 2 z + Hp ( 2 - z) In z] . (B.3) 

Finally we present the coefficient function originating from the Compton 
process (2.17). The asymptotic form is given by (4.12) and can be analytically 
expressed as 
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r4 fl + z 2 



16, 58 > 



-^l„,-f)-2 + 6 2 }ln^ 
3 9 / > iii- 

■l + z 2 



1-z 

m 2 



In 



2 Q 2 + 

— + 11 U ln !-^> 

m z L 1 — z v 3 



+ (t 



2 



-Li 2 1 

3 V 



5 C(2) 



+-ln 2 (l-^) + 41n 2 z- — ln(l 
3 9 

10 



— lnzln(l - 

. 134, 
z) H — — lnz + 



•(2-6z) ln(l -z) + (y - 10s) lnz + 



9 
19 



359 

~27~ 



19^ 



(B.4) 



In the above expression one should bear in mind that the singularity 
at z = 1 will never show up because of the kinematical constraint z < 
Q 2 /(Q 2 +4m 2 ). However after convoluting L^ s '^ by the parton densities, 
the structure function gi(x,Q 2 ,m 2 ) will diverge as ln 3 (Q 2 /m 2 ) in the limit 
Q 2 ^> m 2 . In this limit the upper boundary z max in (2.9) will tend to one 
and the virtual gluon which decays into the heavy quark pair becomes soft. 
The soft gluon annihilation which causes the cubic logarithm above has to be 
added to the two-loop vertex correction containing the heavy quark (Q) loop 
which is calculated in appendix A of [32|. In this way the cubic logarithm 
is then cancelled. The final result will be that in (B.4) the singular terms at 
z — 1 have to be replaced by the distributions (ln fc (l — z)/(l — z)) + defined 
by 



dz ( 



h/fl 



/(*) 



dz 



\n k (l 



{/(*)-/(!)}, (B.5) 



and one has to add the following delta function contribution to (B.4) 



£NS,S+V,(2) /, 



Q 2 



m 



m 



xln( — 



+ fc(2) 



CT / 5(l-,){21n 2 (g)-[f C (2) 



32 



265. 
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Figure Captions 



Fig. 1. One-loop graphs contributing to the OME A$. The solid line indi- 
cates the heavy quark Q. 

Fig. 2. Two-loop graphs contributing to the OME Aq 9 . The solid line in- 
dicates the heavy quark Q. Graphs d.19 and d.20 contain the external 
gluon self-energy with the heavy quark loop. In this loop a sum over 
all heavy quark species indicated by H (m 2 H > m 2 ) is understood. 

PS (2) 

Fig. 3. Two- loop graphs contributing to the OME A Qq ' . The solid line 
represents the heavy quark Q whereas the dashed line stands for the 
light quark q. 

NS (2,) 

Fig. 4. Two-loop graphs contributing to the OME A qq Q . The gluon self- 
energy contains the heavy quark Q with mass m in the quark loop only 
which is indicated by the solid line. The dashed line stands for the 
light quark q. 

Fig. 5. R (5.5) plotted as a function of Q 2 at fixed x; x — 0.1 (upper dotted 
line), x = 0.01 (solid line), x = 10~ 3 (middle dotted line), x = 10~ 4 
(lower dotted line). 

Fig. 6. R (5.5) plotted as a function of x at fixed Q 2 ; Q 2 = 10 (GeV/c) 2 
(lower dotted line), Q 2 = 50 (GeV/c) 2 (middle line), Q 2 = 100 (GeV/c) 2 
(solid line). 

Fig. 7a. g^ xact (Born) (dotted line) and g^ ppiox (Born) (solid line) as a func- 
tion of x for Q 2 = 10 (GeV/c) 2 . 

Fig. 7b. Same as in Fig. 7a but now for Q 2 = 50 (GeV/c) 2 . 

Fig. 7c. Same as in Fig. 7a but now for Q 2 = 100 (GeV/c) 2 . 

Fig. 8a. g[ xact (Born) (dotted line) and ^ pprox (NLO) (solid line) as a func- 
tion of x for Q 2 = 10 (GeV/c) 2 . 

Fig. 8b. Same as in Fig. 8a but now for Q 2 = 50 (GeV/c) 2 . 
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Fig. 8c. Same as in Fig. 8a but now for Q 2 = 100 (GeV/c) 2 . 

Fig. 9a. g\ ght (solid line) and g^ pprox (dotted line) both in NLO as a function 

of x for Q 2 = 10 (GeV/c) 2 . For x < lO" 3 both #} ight and #i Pprox become 
negative so that we have taken their absolute values. 

Fig. 9b. Same as in Fig. 9a but now for Q 2 = 50 (GeV/c) 2 . 

Fig. 9c. Same as in Fig. 9a but now for Q 2 = 100 (GeV/c) 2 . 
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Fig. 2 ( continued ) 




Fig. 4 
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Fig.7a 




Fig.7c 




Fig.8a 




Fig. 8b 




Fig.8c 




Fig. 9b 




Fig.9c 



